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Abstract 

We determine the UHE neutrino-nucleon cross-sections analyti- 
cally both for charged and neutral current in the leading order by 
using the solutions of non-singlet DGLAP equation as well as singlet 
DGLAP equation respectively. Next we determine the cross-sections 
for both charged and neutral current at one loop level by earring 
out numerical integration of the double differential cross-section. For 
analytical solution, we take the standard MRST 2004 f4 LO parton 
distributions at Qq = 1 GeV^ as input and for numerical solution 
(exact), we consider the GRV 94 parton distributions. 

Our analytical as well as numerically determined results are com- 
pared with the NLO results of various authors. We find our numer- 
ically determined result in sufficient agreement with the results ob- 
tained by those authors. On the other hand, our analytical result 
matches better with other results at lower end of the energy spectrum 
of the UHE neutrino rather than the higher end. 
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1 Introduction 



Charged current induced interactions in the Ultra-High Energy (UHE) neu- 
trino (antineutrino)-nucleon scattering processes are represented as 

v^N — !■ /i~X, V^N — )■ fi'^X, where = {n + p)/2 is an isoscaler nucleon. 

On the other hand, neutral current interactions are represented as: 

In case of charged current interactions, the mediating particles are W- 
bosons (W^^) and in case of neutral current interactions, the mediating par- 
ticles are Z-bosons {Z'^). 

Study of UHE energy neutrino-nucleon interaction on the basis of stan- 
dard model of particle physics has been pursued by various authors [H [2l [3l IH 
El El [7] from time to time since long back. Recent measurements indicate that 
there are GZK-violating air showers with reliably determined energies above 
10^^ GeV [HI E] triggered off by UHE neutrinos. Event rates of these showers 
can shoot up the UHE cross-sections several thousands times as predicted 
by standard model. This unusual bumping up of cross-sections for neutrinos 
having energies above 10^ GeV gives rise to a puzzle quite unexplainable by 
the standard model. 

UHE neutrinos having energies above 10^ GeV collide violently with nu- 
cleons at center-of-mass energies [10] which is greater than 



At such gigantic energies, the magnitudes of the ultrasmall Bjorken x 



= V^MK - 14 GeV TeV 
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values are less than 

' \ /0^\ / 10« GeV \ 

M^J [y )[ I 

At such high energies several approaches [101 [H] based on physics beyond 



X ~ 2 X 10^ 



the standard model have been proposed to explain the event rates measured 
by neutrino telescopes hke 'Pierre Auger Observatory' [121 [IHl Ell [HI 116] . 
ICECUBE [I711T811T9], 'OWL' [20j, 'EUSO' [21] and 'JEMEUSO' [22]. One 
such approach is based on electroweak instanton-induced processes where 
there is interaction between cosmogenic neutrinos and the nucleons in the 
atmosphere. Such a process, which was proposed by Ringwald [23] and 
Bezrukov et al [211 125] violates the baryon+lepton number, B + L. Other 
approaches include production of miniature black holes in ultra-relativistic 
collisions with energies above the fundamental Planck scale Mr, [Md ~ TeV) 
[251 l26l [271 [281 [29] , neutrino interactions through TeV string resonance pro- 
cesses with string scale Ms [30j, production of quasi-stable charged sleptons 
(staus) predicted by weak scale supersymmetric models with the supersym- 
metry breaking scale greater than 5 x 10^ GeV [311 [321 [33] . Stau, which has 
a decay length on the scale of 10 km, is the next-to-lightest particle (NLP) 
which is produced by neutrino interactions on earth or atmosphere. Detec- 
tion of pairs of charged tracks of secondary particles produced by staus in 
neutrino telescopes like 'ANITA' [M] [35] would be an excellent way to probe 
the SUSY breaking scale. A relative comparison of neutrino-nucleon inter- 
action cross-section as per standard model and beyond is observed in Figure 

m 
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QCD I 
EW instaiiton 

QCD with saturation "| 
black hole (M=lTeV, M'"°=5TeV, n=4) : 

le+07 le+08 le+09 le+10 le+11 le+12 

E[GeV] 

Figure 1: Neutrino-nucleon interaction cross-section as per standard model 
and beyond. Figure taken from Ref. [11] 

In the unpolarized case, several authors had solved DGLAP evolution 
equations [361 133 EH EH HQ] efficiently even up to NNLO accuracy with 
numerical technique. There are basically two approaches: one by using the 
Mellin tranform method (N space based) and the other by using the x space 
method. The numerical efficiency of the ffist method is more than the sec- 
ond, but the second method is more flexible than the ffist because here one 
requires inputs only in x space. Now-a-days, one can get several numerical 
packages publicly to solve the DGLAP evolution equations efficiently and ac- 
curately, e.g., PEGASUS [H] , HOPPET QCDNUM Pj, CANDIA 
[i5| 116] . For example, one can use the Fortran package QCD-PEGASUS 
to perform the evolution using the symbolic moment-space solutions on a 
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one-fits-all Mellin inversion contour to get a flexible and very accurate solu- 
tions of the evolution equations. PEGASUS implies Parton Evolution Gen- 
erated Applying Symbolic U-matrix Solutions'. On the other hand, there are 
several high flexible x-space methods. For example, the HOPPET NNLO 
parton evolution package is a very good Fortran package for carrying out 
DGLAP evolution and other common manipulations of PDFs. Similarly, 
the evolution equations for parton densities and fragmentation functions in 
perturbative QCD can have numerical solution with the help of the QCD- 
NUM program, which is powerful to evove un-polarised parton densities upto 
NNLO in powers of the strong coupling constant as- Precision Studies of the 
NNLO DGLAP Evolution at the LHC can be carried out with programme 
CANDIA (progrmming language C and Fortran) [45j. Similarly, NNLO Log- 
arithmic Expansions and Precise Determinations of the Neutral Currents 
near the Z Resonance at the LHC (the Drell-Yan case) can be obtained us- 
ing CANDIA [36] . In spite of having these highly efficient numerical packages 
of solving DGLAP equations for the QCD evolution of parton distributions 
(PDFs), analytical approahes for solving DGLAP equations bear special in- 
terest. It is interesting to see whether the analytical solutions deviate at all 
from the exact numerical solutions and if so, to what extent. 

In recent years, an approximate method [371 I3HI HHl [50l [51] to solve 
DGLAP equations [361 ISH (Ml [391 130] has been pursued with considerable 
phenomenological success. In that approach, we have used Taylor series 
expansion to express these equations as partial differential equations in x 
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X = 



2p.q 



) 



and t [t = log — , valid at low x. 




The aim of the present paper is to apply these approximate solutions at 
low X to the Standard Model Interactions of ultra-high energy neutrino with 
nucleon which has an important bearing in neutrino astronomy. In doing so, 
we assume that the approximate solutions of DGLAP equations are valid in 
the region of ultrasmall values of Bjorken x. With this assumption, we will 
find leading order (LO) expressions for UHE cross-sections for charged as well 
as neutral neutrino-nucleon interaction [ctcc, o'nc) within this formalism. All 
throughout our earlier papers [521 [53] we have resorted to MRST distributions 
and obtained results on the basis of these distributions in non-singlet as well 
as singlet case. We are going to apply these results in the present paper, so 
here too we consider MRST distributions [54l |55] to determine our analytic 
results. 

We will also find UHE neutrino-nucleon cross-sections by using the nu- 
merical solutions of DGLAP equations. For this purpose, we take the help 
of GRV94 parton distributions [56] . 

There are several authors [H [21 El [5l [571 [HE] who has predicted the UHE 
neutrino-nucleon cross-sections by different techniques. We will then compare 
our predictions with the result of several such authors and study the range 
of validity. In Section [21 we outline the formalism while Section [3] contains 
results and discussions and Section [D contains comments and conclusions. 
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2 Formalism 

2.1 Approximate solutions of DGLAP equations: 

The approximate solutions of DGLAP equations |36l [STJ [381 EHl HQ] for non- 
singlet structure function F^^{x,t) and singlet structure function F2{x,t) 
are respectively given by 



/ f \ n{x,t) 

t) = to) i^-J (n(l, t) > 0, t > to 



(1) 



and 



Fi{x,t) = F'{x,to)^ 



X^ix] 



k(x,t) 



-\k{x,to) 



(A;(l,t)>0,t>to) (2) 



where t = log (^^^ and to = log(^^^, Q^(= — q'^) is the four momentum 
transfer and is the QCD cut-off parameter. X^{x) is given by |52] 

dx 



where 



X (x) = exp 



Pi{x,t) = -AfX 





f (1-1=) 







(3) 



(4) 



It is to be noted that in DIS neutrino scattering, F'^^{x,t) is identified as 
xF3{x,t). 
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2.2 UHE neutrino-nucleon cross-section 
2.2.1 Charged Current- Analytical Solution: 

The differential cross-section for the charged current (CC) UHE interaction 
is given by [3], 



dxdy 2tt \ 



(5) 



where Fermi coupling constant Gp = 1.1663 x 10~^GeV~'^, s = 2ME^(y^ 
is the centre-of-mass energy squared, x is the usual Bjorken variable which 
represents fraction of the nucleon's momentum carried by the quark and 
y = — , where u = E,y — = Ehad- Ey represents the energy carried by 
incoming neutrino, E^ represents the energy carried by outgoing muon and 
Ehad represents the energy carried by the final state. 

The total cross-section for UHE CC interaction in an isoscaler target 



A^[= -^^i^] is given by. 



^1 ^1 j2^^{v)N 



Replacing y by in the L.H.S. of eq. |5l we get the expression for the 
double differential cross-section for charged current as 



cc 



Q2 



dxdQ'^ 




X, Q') ± Y^xF^ 



x,Q' 



(7) 
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where 



(8) 



and 



?CC _ T7CC 



2xFl 



cc 



(9) 



The total cross-section for UHE CC interaction can be re- written as [3 [59] 



a. 



u{y)N 
CC 



dx / dQ 







'CC 

dxdQ"^ 



(10) 



At one loop level (LO), the Callan Gross Relation Fi = {^)F2 holds. So 
obviously, Ff^ = 0. 

Hence, putting the value of 1± in eq. ([7]) and considering t = log (^^^ as 



in Subsection I2.H we get 



d'a 



CC 



dxdQ^ 



27 



1 + 



M^j X 



Using eq. ([T]) and eq. (I2D, eq. ffTT]) can be re- written as. 



(11) 



d'a. 



CC 



dxdQ^ 



2n \ 



k(x,t) 



1 

X 



l-y + ^]F^^''^ 



CC f t\ 



±y{l-'-)xF', 



where k{x,t) is identical with H^^p^{x,t) given by 



(x, tn) I — 

tr 



^ \ n(x,t) 



(12) 



H^{x, t) = 1.284[-4.64x + 1.02(1 - a;)](1.52 - 1.16s - 0.06s^ 



(13) 
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and n{x, t) is identical with H{x, t) given by [52] 

H{x, t) = 0.380[-5.796x + 0.996(1 - x)] (14) 
The v{v)N structure function Fi can be written as; 

= fI'"''' + Fl^ (15) 

where i^^'*^'** provides contributions of all light (u, d, s) quarks and F^ 
refers to contributions of all heavy (c, 6, t) quarks. 

We have the following LO expressions in QCD for charged current (CC) 
interactions [3], 

pu,Ught _ 2^jpi',light ^^^^ 
pu,Hght ^ _ + 3) + + ^)|^^^|2 ^ 2s|K.P (18) 

where m = u{x, Q^) etc. and Vud, Vus etc. are the relevant CKM matrix 
elements [60| WTj 162] . We have the following values for the following matrix 
elements: Kd = 0.97419 ± 0.00022, Ks = 0.2257 ± 0.0010. Here we consider 
the individual contributions from the CC subprocesses of the light quarks 
i.e. W + d^u, W + s^u, W + u^d, etc. Let us now consider the 
heavy quark treatment in a fully massive way [3] as we find in case of GRV 
98 parton densities [63j and have the following expressions for heavy quark 
component. 
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puMavy ^ ^ ^^|^^^|2 ^ 2xs\V,s? + 2xC + ^ ^2)|^^^|2 (20) 

p-Mavy ^ (^ ^ ^)|^^^|2 ^ 2s|F,,|2 _ 2c + 26(e, + (21) 



2 ^ 

mf 



Here ^ = x ( 1 + -|- , = 175 GeV, nib = 4.5 GeV, = 0.2256 ± 



2 



0.0010 [62] and contributions like — |- have been neglected. The Charm con- 



tribution to the heavy flavour component F^'^'^'"^ at one loop level takes place 
through the transition channels W^s — )■ c and W^d — )■ c and the top bot- 
tom contribution comes through the fusion subprocess channel W^g — i- tb. 
There are also additional subprocesses W'^g — )■ cs and + s' — )■ gc, etc., 
where s'^j^ = \ Vcs\'^s+ \ Vcd\'^{d + u)/2 [M]. Again we neglect the contribution 
of t{x, Q"^) since at the relevant range in the UHE domain, the nucleon 
contains a very negligible amount of ti sea [57]. It has been found that the 
contributions from the cs sector surpasses the minor contributions provided 
by heavy tb in LO where we take the fully massive m^^t ^ 0. This situation 
prevails for E^, < 10^ GeV, but when we progressively move towards higher 
energy and make a switchover from LO to NLO, then we find that LO con- 
tribution gets an additional 0{as) correction by about 20% due to the NLO 
massive tb contribution ^4J as evident in Figure |2l 

On the other hand, if we consider heavy quark flavours as massless, as 
we find in case of GRV 92 parton densities [65] 166] 



p.Mavy ^ 2c(v;^ + Al) + 2b{Vj + Al) (22) 

pu,hea.yy _ 2j.p^^^'^°-'"y (23) 
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10- 10^ lO"* ICr' 10*' 10^ 10^ lO' 10'° 10^' 10 

E„(GeV) 



Figure 2: Total Neutrino- nucleoli interaction CC cross-section versus Energy 
of the neutrino at NLO along with individual contributions due to light, 
charm quark transitions as well as contributions due to fully massive tb. 
Figure taken from Ref. [HI] 

pu,heavy ^ q ^24) 

where K = | — | sin^ 9w, = —^d = | and sin^ 9w = 0.232. 

Now considering the minimal contributions provided by the heavy quarks, 
which is not so prominent as the contributions provided by the light quarks, 
we only focus on the light quarks in our work. 

For uN interaction, one uses the relations eq. f|T6|) to eq. f lTSj) and the 
following correspondances, 

FI2 = Fl^{q ^ q) (25) 
F^=-F}^{q^q) (26) 
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Now putting the explicit parton distributions from eq. ( 1TB]1 to eq. ( TTSi) 
in eq. (fT2!) . we have; 



cc 



dxdQ"^ 



2^ 



1 + 



X 



l-y + yj2a;|^ \u{x, to) + d{x, to) 



±2/ (^1 - 1^ X {- to) + d{x, to) 

+ [dix,to)+uix,to)]Kd\^ + 2s\Vus\^} i^-J (27) 

From the kinematics of deep inelastic neutrino scattering, we have 

S-M2 

E,. -- 



k{x,t) 



2M 



(28) 



In the UHE region of interest scanned by us, we can reasonably re-write eq. 
(EHl) as [a [58] 



(29) 



2M 2M 

Again from the kinematics of deep inelastic neutrino scattering, we can 
reasonably approximate O [58] 



y 



(30) 



x{s — M^) xs 

By setting proper lower cuts for and x integration [5l [58] in eq. ([10 

we get 



a. 



v{v)N 
CC 



dx / dQ 



91 Ql 



2 " '^cc 
dxdQ^ 



(31) 



13 



The contribution from xF^ at Ultra- low Bjorken x in UHE domain is very 
small. But still retaining this factor, we have from eqs. ( 127|) and (I3T]) . 

-2 1 



0", 



u{V)N 

cc 



x2x <^ - 
.2 



w . 



dx 

X 



1 - 



xs Ix^s^ 



u{x,to) + d{x,to) + - [d{x,to) + u{x,to)]\Vud\^ + 



X^ix) 



k{x,t) 



t, 



k{x,to) 




2xs' 



)x |- u{x,to) + d{x,to) 



■) / t \ "(^'*) 
+ [d{x,to) + uix,to)]\Vj^ + 2s\Vus\''\ l^- 



(32) 



Here u{x,to), d{x,to) etc. represent parton distributions at various values of 
X and input value to- 

2.2.2 Charged Current-Numerical Solution (exact): 

Putting the explicit parton distributions from eq. flT6|) to eq. f|T8|) in eq. 
(ITT]) , we get the expression of double differential cross-section as 



d^a. 



cc 



dxdQ'^ 



2^ 



Ml;, 



-2 



t) + d(x,t)]) 

+ ^[d(a;, t) + M(x,t)]|Kdr + s|K.r} ± y (^1 - I) X 
{- [M(a;,t) + d(x,t)] + [d(x,t) +M(x,t)]|Kdr + 2s|Ksr} 

(33) 



Now putting the expression for double differential cross-section from eq. 
in eq. ([SI]), we get 
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i/(F)iV 

cc 



Ql ^ 



-2 1 



dx 

X 



+d{x,t)] + - [d{x,t) +u{x,t)]\Vud\'^ + s\Vu 



,21^0! 

XS 



2xs 



X 



{-[u{x,t) + d{x,t)] +[d{x,t)+u{x,t)]\Vud\^ + 2s\Vus\^}] (34) 



Here u{x,t), d{x,t) etc. represent parton distributions at various values of x 
and t. 

2.2.3 Neutral current- Analytical Solution: 

The total cross-section for UHE NC interaction {aNc) in an isoscalar target 
A^[= is similar to eq. (^T]) . viz 



where 



u(u)N 



dx / dQ 



2U_Ojfc_ 

dxdQ"^ 



(35) 



,2 v(F)7V 



NC 



dxdQ"^ 



G 



271 



Ml;) X 



1 -2/ + IT 1^2 



(a;, to) ( -j 



X 



fc{a;,t) 



(36) 



[xs{x)r 

to be compared with eq. (fT2|) for charged current. 

We have the following LO expressions for structure functions involving 
Neutral currents (NC) [3]: 
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(37) 

^u,iight ^ 2xFr'''^^* (38) 
^p..,iight ^ 2(^^ + 4)(KA« + VM (39) 

where the expressions for Vu and Vd as well as the values of Au and Ad have 
been mentioned beforehand. 

As in the case of charged current, we focus only on the contributions from 
light quarks. 

Putting the values of V^, A^, Vd and A^ in eqs. fl371) . fl38l) and fl39|) and 
using eq. (1361) and eq. (l35l) . we finally get; 



(T 



u{u)N 
NC 



27r 



rfgM 1 + 



Ml, 



Q2 



xs 2x^s^ 



{0.1432X 



X 



+ 0.1849X 



m(x, to) + ^t(a^, ^) + d{x, to) + (i(x, to) 
u{x, to) + u{x, to) + d{x, to) + d{x, to) + 4s | 



k(x,t) 



xs \ 2xs 



^ \ n(x,t)' 



tc 



|o.268x u{x,to) — u{x,to) + d{x,to) — d{x,to) | 

(40) 



which can be re-expressed as 
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u(V)N 
NC 



G 



F 

27r 



-2 1 



Ml 



dx 

X 



xs /x^s^ / 



s 

X [{uy{x, to) + dy{x, to)) + 2 {u{x, to) + v{x, to))] + 0.1849x 
X [{u,,{x, to) + dy{x, to)) + 2 to) + v{x, to) + 4s)]} -j^^ 

tn 



X 



X^{x) 



k(x,t) 



{0.268x[uy{x,to) + dy{x,to)]} [j- 



2.2.4 Neutral Current-Numerical Solution (exact): 



(41) 



We have in case of neutral current the following equation 



dxdQ"^ 



Ml, 



27r 



±y ( 1 - 1 ) xf; 



X 



{x,t) 



x,t) 



(42) 



similar to eq. ffTTl) for charged current. 

Now putting the explicit parton distributions from eq. (157|) to eq. (15^ 
in eq. ( H2|) and then using eq. (135|) . we get the expression for neutral current 
cross-section as 



a 



u{u)N 
NC 



Ml 



-2 1 



Q2 



dx 

X 



xs 



2x^s^ 



{0.1432a; 



X 



X 



u{x, t) + u{x, t) + d{x, t) + d{x, t) + 0.1849X 



u{x, t) + u{x, t) + d{x, t) + d{x, t) + 4s 
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} 



± 



xs 



2xs 



{0.268x [u{x, t) - u{x, t) + d{x, t) - d{x, t)] } 

(43) 



which can be re-expressed as 



2 1 



u{v)N Gp f 2 / 1 , \ f dx 



xs ~^ 2x^s^ j ^ 



X [Mx, t) + dy{x, t)) + 2 t) + v{x, t))] + 0.1849X 



X t) + d.,{x, t)) + 2 (Il(a;, t) + v{x, t) + 4s)]} 

xs \ 2xs , 



±9_ h _ ^1 {0.268a; t) + t)]} 



(44) 



2.3 Role of the gauge-boson propagator: 

For neutrino energy E^, <^ -^-^ (or say Q"^ ^ ^w^)? it is found that a (x E^. 
This is because of the fact that in the hnear regime at fow Q^, the gauge boson 
propagator and the PDFs do not almost depend on and then for Eyyi^f' 
GeV, 0" is found to be approximately proportional to When E^ is greater 
than 10^ GeV i.e. for E^ > (or say Q'^ > M'^), the gauge boson 
propagator decreases sharply making the double differential cross-section to 
die off fast [2]. Consequently the gauge boson propagator puts a restriction 
on to assume values near and thus integration region shrinks, thus 
making the upper limit of the integral proportional to M^r O . Thus 
at ultra high energy, dependence of the gauge boson propagator and the 
PDFs becomes prominent. In one hand, there is a power suppression from the 
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boson propagator, while on the other hand there is a logarithmic growth of 
the PDFs. In overall, the propagator dominates and it leads to the generation 
of the dampness of the total cross-section [68J. It is pertinent to note that 
the effective interval in the fractional parton momentum x gets limited to 
the region around . It is interesting to note that in the (x, Q^) phase 
space at the high energy of the tone of E^, > 10^ GeV, the cross section 
gets dominated by the behavior of the quark densities at x = ^ and 

~ [69]. For the neutrino energy E^^ ~ lO"*^^ GeV, we have x ~ 10 ^. 

Though a cut off in Q"^ is automatically created at Q"^ ~ M^^ due to the 
propagator effect, but we have carried out the integration upto the upper 
limit s = 2MEiy (eq. (1291) ) to take into consideration the unaccounted part 
of the results of cross-section for > M"^. Of course, we have once checked 
that the difference is not too much, rather minimal by putting the upper 
limit at ~ ~ 10'* GeV^, thus justifying the argument of the validity 
of the taking the upper limit of the Q-integration around ~ at ultra 
high energy of the neutrino (as well as antineutrino). 

2.4 Gluon recombination, saturation and unitarity: 

It is interesting to note that in these approaches, cross-section increases 
sharply with Ey in a fast power-law-like growth at the highest energies. But 
this rapid growth gets slowed down when we consider the upper extreme end 
of the low X region and unitarization of the cross section takes place [H [TJ |67]. 
At extremely low values of x, the density of parton (gluon) phase space will 
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-1 - 



saturation region 




Figure 3: The figure in the above is the contour plot of the neutrino-nucleon 
cross section ( ^^2*^^ ) in the (x, Q^) plane. Saturation region is plotted from 
Eq. psjl . Figure taken from Ref. [67] 

increase. Consequently, there will be overlapping in transverse space and 
gg ^ g recombination processes will trigger up. The characterizing scale of 
the saturation region (regime in x and Q"^) is given by, 



The value of the parameters A and xq are obtained as 0.288 and 3.04 x 10 
respectively by a fit to HERA data. The domain of saturation starts with 




(45) 
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those values of Q which are lower than the saturation scale Qs{x) at a given x. 
It is clear from the above equation that with the increasing/decreasing values 
of x, saturation takes place at smaller/larger values of Q^. In the highest 
sensitive UHE neutrino interaction domain x ~ 10~®, saturation scale is 
around Ql = 20 GeV^. It is evident from the Figure [3] that total cross-section 
is dominated by the scale Q ~ Myyz- An obvious question is that what is 
the effect of gluon recombination on the total cross-section? Interestingly, 
the contribution of the saturation region to the total cross-section is not so 
much appreciable [U EZ]- Even at very high energy say E^, = 10^^ GeV, 
the contribution is much lower than 1%. The region of < 1 GeV^ and 
consequently the saturation region at E = 10^^ GeV contribute very little 
to the total cross-section. We must realize that at the time of evaluation of 
total cross-section, we carry out our integration over large range of values of 
X and in the (x, Q^) phase space and since the saturation region comprises 
of a very small section of the phase space, hence the effect is marginal. But 
the importance of saturation lies in retarding the growth of cross-section so 
that it does not cross the unitarity bound (Froissart bound). This is quite 
obvious from the Figure |3l 

2.5 Small x extrapolation 

Since the UHE neutrinos coming from various sources have in general en- 
ergy more than 10^ GeV, so perturbative calculations should be carried 
out using parton distribution functions at Ultra- low Bjorken x {x < 10"^). 
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Particularly, the sensitivity of UHE neutrino interactions are greatest in 
the domain x ~ 10^^ and ~ M^r. But high energy measurements 
of deep inelastic scattering (DIS) of lepton-nucleon at the collider HERA 
UniiniESlESlCllESlESlEIlEH] at DESY have provided interesting data 
from time to time. To the best of its ability, it has provided data for x > 10^^. 
The greatest available energy ^/s^ available at HERA is 319 GeV, which 
is far below the expected UHE neutrino-nucleon collision energies of about 
y/SuN = 10*' GeV [31 . So it is not possible to scan the ultra-high energy 
region and evaluate the UHE cross-section from the present DIS experiments 
alone. 

There are two different approaches to determine the cross-sections for 
UHE neutrino-nucleon interaction. One way is to parametrize the parton 
distribution functions by global fitting to a rich set of available data at low 
values of (say = Ql). The numerical solutions of the DGLAP equa- 
tions [361 E3 EHl EH SO] are then used to evolve parton densities at higher 
values of (Q^ > Ql)- This Q^-evolution of small-x parton distribution is 
not the right choice to scan the region of UHE neutrino-nucleon interactions 
as the evaluation is valid in the region of x > 10^^. Extensive extrapola- 
tion [21 [31 m EH EHl Eni Ell ES] is then carried out to determine the parton 
densities from the (x, ) domain of terrestrial accelerator HERA to the 
ultrasmall values of Bjorken x, where there is no data. But unfortunately, 
this type of extrapolation brings about some uncertainties in the estimated 
UHE neutrino-nucleon cross-sections [3] measured with the help of neutrino 
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telescopes. This is because we are not so sure about what is happening in the 
region x — )■ and our ignorance forces us to take different assumtions ( or 
educated guesses) about the behaviour of the distribution functions at ultra- 
small x, leading to the large variations of different evaluated cross-sections. 

As an example of the first approach, we can mention the case of Gandhi et 
al P EZl EH ED] who used CTEQ3 and CTEQ4 parton distributions [SSI El] 
and MRS parton distrbutions [851 ESI EH [88] ( MRS A',G,D.^ and D'_) 
to determine the charged current cross-section acc and the neutral current 
cross-section a^c- On the other hand, in the second approach, Gluck, Kretzer 
and Reya [3] applied QCD- inspired radiative parton model [561 ESI ESI EH EH 
ED] to evaluate parton densities directly for x < 10^^. The model utilizes the 
principle of QCD dynamics to extrapolate parton densities to kinematical 
regions as yet unexplored by the present day terrestrial experiment and to 
go for respective fits. Parton densities obtained within the framework of this 
model only are claimed to be free from the ambiguity of the earlier approach. 
Here one calculates the parton densities in the region x — t- by the direct 
application of QCD dynamics, without introducing any fit parameter in that 
region. It may be mentioned that though MRS and CTEQ parametrisations 
do not work at all for values of x below 10"^, but GRV parametrisation works 
well for values of x upto 10"^ p. 

At Ultra-low Bjorken quarks will hugely multiply in a dominating 

way and these arise from gluon splitting into quark- ant iquark pairs. So the 
technique of gluon extrapolation will provide one a direction to carry out sea 
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quark extrapolation. If the gluon extrapolation is carried out at some input 
value Qq in the following way [M] . 

xg{x,Ql)^ AiQDx-"- a; « 1 , (46) 
then at higher values of Q^, we can extrapolate gluons similarly 

xgix.Q'') ^ A{Q^)x-^ . (47) 
We can extrapolate sea quarks in a similar approach 

a;g(x,Q^) = ^-^^ a;g(a;™„,Q^) . (48) 

One can determine A from the PDFs at = for each flavor. It is 
pertinent to know that we have developed our formalism in this paper in 
leading order only, so we will use LO parton densities in this case. In case 
of NLO, one has to use NLO parton densities strictly, with the appropriate 
addition of convolutions with the fermionic Wilson coefficient Cq and the 
NLO contribution Cg® g [93j, though these additional contributions never 
provide more than 2% contributions [3]. As a simplified case, we have not 
gone for small x extrapolation in this paper with the assumption that there 
is not too much remarkable difference between the extrapolated and the non- 
extrapolated results. 

2.6 Equivalent neutrino-nucleon cross-section 

The two detector collaborations HI and ZEUS at HERA at DESY in Ham- 
burg recorded integrated luminosity of 1 ph~^ from May to October, 1993. 
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Figure 4: Neutral and charged current cross sections as measured by HI and 
ZEUS alongwith Standard Model predictions. Figure taken from Ref. [M] 

The amount was twenty times the luminosity that was recorded in 1992 dur- 
ing its run for the first time, where the energy of the electron beam was 26.7 
GeV [95j. In 1993, energy of the electrons was raised to 27.5 GeV, which were 
allowed to collide with 820 GeV protons, where the centre of mass energy 
(v^) was around 300 GeV. From July 18, 1994, electron beam was replaced 
by positron beam due to the limitation imposed by the lifetime of the elec- 
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Figure 5: Equivalent charged current cross sections in deep inelastic neutrino- 
nucleon scattering as measured by HI. Figure taken from Ref. [95]. Similar 
figure is also obtained in |96] 
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tron beam [M]. In the HERA II phase in 2000-2001, the collider as well 
as experiments were upgraded and consequently 920 GeV proton beam was 
allowed to collide with 27.6 GeV electron/positron beam with the centre-of- 
mass energy 319 GeV and integrated luminosity around 150 pb~^ [51]. Recent 
results from HERA shows HI luminosity as 181 pb~^ in e~p and 294 pb~^ in 
e~^p scatterings respectively [57] . 

On the basis of data collected by HI and ZEUS, the neutral and charged 
current cross-sections were calculated during the HERA I phase (1994-2000) 
and was found to coincide nicely with the electroweak Standard Model predic- 
tions based on CTEQ6D parton density functions [91] as revealed in Figure 

SI 

From neutrino-nucleon scattering experiments in the eighties, weak charged 
currents were studied, where it was found that the cross-section is directly 
proportional to neutrino energy [98]. This linear behaviour of cross-section 
clearly pinpoints the fact that the sensivity of the neutrino energy to the 
W Propagator effect is extremely small and hence negligible. The field of 
electroweak physics was opened for the first time when charged current re- 
action ep — )■ uX was studied at HERA, where the role of W Propagator 
became important becauge of the high momentum transfer comparable to 
Q"^ = = {SOGeVy and the equivalent fixed target energy of the HERA 
collider was 50 TeV. Since this process is just the inverse of neutrino-nucleon 
scattering, the measured ep cross-section could easily be converted into an 
equivalent uN cross-section at 50 TeV [95l ES] . Interestingly, this value of 
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equivalent cross-section comes very close to the value of cross-section as pre- 
dicted by Gandhi et al in the recent past. In Figure [5l low energy neutrino 
data has been represented by the crosses, whereas the full and open stars 
represents the equivalent vN cross-section as obtained from the HERA ex- 
periment. The extrapolation from low energies has produced the straight 
line at higher and higher energies assuming M^y = oo, but when the effect 
of W propagator is considered, one gets the bent curve which represents the 
predicted cross-section. The vN cross-section corresponding to equivalent 
fixed target energy of 50 TeV is represented by HI point in the Figure O 

3 Results and discussions: 

To estimate the cross-sections for charged and neutral current from eqs. (1311) 
and we use MRST 2004 f4 LO input parton distributions at = ii\q = 
1 GeV^ [511 [55] 

Putting these distributions in eqs. ( 132|) and ( HTj) and using eqs. ( l3Ti) and 
( !35|) to integrate over x and Q^, we obtain the values of cross-sections for uN 
(as recorded in Table [1]) and uN (as recorded in Table [2]) interactions both 
for charged current and neutral current respectively . 

As expected, for 10 < < 10^ GeV, the cross-section rises linearly with 
energy as discussed in Subsection 12.31 and 12.61 

Then the cross-section versus energy relation have the following power 
law forms: 
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(J, 



vN 
CC 



a 



uN 
NC 



1.70 X 10-35 
2.24 X 10-35 

3.18 X 10-36 
5.28 X 10-36 cm^ 



^^g^. 0.199 

1 GeV) 
I GeV J 



0.183 



vN 
CC 



1 GeV J 
Eu \ 

I GeV J 

1.02 X 10-35 cm^[ ^ \ 
GeV J 



0.213 



0.185 



0.231 



2.23 X 10-35 cm2 



Et7 \ 

I GeV J 



a 



VN 
NC 



1.97 X 10-36 cm 



2 1 ^ \ 
1 GeV) 

5.28 X 10-36 cm^[ j-^^) 



0.183 



0.239 



0.185 



10^ <E^< 10^ GeV 
10^ <E^< 10^2 Gel/ 

10^ <E^< 10^ GeV 
10^ <E^< 10^2 Gel/ 

10^ <Ejj< 10^ GeV 
10^ < < 10^2 GeV 

10^ < < 10* Gel/ 
10* <Etj< 10^2 Gel/ 



(49) 



(50) 



(51) 



(52) 



Next we use GRV94 parton distributions [56j for determining the cross- 
section numerically. We use Mathematica 4.1 version to carry out numerical 
integration with appropriate limits. To parametrize LO parton distributions 
as per radiative (dynamical) predictions, one has to define 



S = l0£ 



log 


[ 1 


(0.232 GeVY_ 


log 


Ac 


(0.232 GeVY_ 



(53) 



to be evaluated for = 0.23 GeV^. The validity of all the following 
parametrizations are in the range 0.4 < < 10^ GeV^ (i.e. 0.3 < s < 2.4) 
and 10"^ < X < I. 
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The non-singlet structure functions can be parametrized as 

xv{x, = Arx"(l + Ax^ + Bx + Cxi)(l - x)^. 
For V = m,. 



a = 0.590 - 0.024s, h = 0.131 + 0.063s, 
N = 2.284 + 0.802s + 0.055s^ 
A = -0.449 - 0.138s - 0.076s^ 
B = 0.213 + 2.669s - 0.728s^ 
C = 8.854 - 9.135s + 1.979s^ 
D = 2.997 + 0.753s - 0.076s^. 

for V = d„ 



a = 0.376, b = 0.486 + 0.062s, 
N = 0.371 + 0.083s + 0.039s^ 
A = -0.509 + 3.310s - 1.248s2, 
B = 12.41 - 10.52s + 2.267s^ 
C = 6.373 - 6.208s + 1.418s^ 
D = 3.691 + 0.799s - 0.071s^ 



The gluon and sea u + d distributions are parametrised as 



xw{x, Q'^) 



X 



(A + Bx + Cx^) + (\og - 
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X exp -E + \ E'sl^ log- 



X 



[l-x) 



D 



where iox w = u + d 



a = 1.451, ^ = 0.271, 
a = 0.410 - 0.232s, b = 0.534 - 0.457s, 
A = 0.890 - 0.140s, B = -0.981, 
C = 0.320 + 0.683s, 
D = 4.752 + 1.164s + 0.286s^ 
E = 4.119 + 1.713s, E' = 0.682 + 2.978s. 

Parametrization of the strange sea distribution is as follows: 



s 

xw'{x,Q^) = -y, (l + Ay/^ + Bx) {1 - xf 



X exp I + W E'sl^ log - 



i.e for ^u;' = s = s 



a = 0.914, /3 = 0.577, 
a = 1.798 - 0.596s, 
A = -5.548 + 3.669x/s - 0.616s, 
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B = 18.92 - 16.73^5 + 5.168s, 



D = 6.379 - 0.350s + 0.142s^ 



E = 3.981 + 1.638s, E' = 6.402. 



(60) 



Putting these distributions in eq. f l34j) and eq. ( I44j) . we obtain the values 
of cross-sections for uN (as recorded in Table |3]) and VN (as recorded in Table 
H]) interactions both for charged current and neutral current respectively. As 
in our analytical case, the following power law forms indicate the cross-section 
versus energy relation: 



/ jt; \ 0-433 
1.05 X 10-36 cm2( — ^) 
uN _ t \l GeV J 

3xl0-c^^(^) 



10^ < E,, < 10^ GeV 



10« < E„ < GeV 



(61) 



cr 



/ TP \ 0.426 

2.34 X 10-37 cm2(— ^) 
_ J VI GeV J 

Nc — S / E \ °'3*'7 

5.68 X 10-37 cm^l — ^ ) 
VI GeK / 



10^ < E,, < 10^ GeV 



10« < E,, < 10^' GeV 



(62) 



(J, 



0.466 



/ E— \ 

6.05 X 10-37 cm^l 

VN _ ) \l GeV 

3xl0-36c.n^(^) 



10^ <E^< 10« GeV 



10« < < 10^' Gel/ 



(63) 



NC 



7.90 X 10-38 cm2 



5.68 X 10-37 cm^ 



0.487 



1 GeV 

_E^y- 

I GeV J 



10^ < E^< 10^ Gel/ 



10« < ^77 < 10^' GeV 



(64) 



This is to be compared with NLO expressions for cross-sections of Ref. 
[3], calculated on the basis of GRV 98 parton distribution [63] 
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vN 

cc 



1.10 X 10~3^ cm 



^ 0.454 

,1 Gev) 

5.20 X 10-36 cm^f — ^ ) 
,1 Ge\/ / 



10^ <E,< 10^ Ge\/ 



10« < < 10^' GeV 



(65) 



/ \ 0-467 

, 3.55 X 10-37 cm^—f-) 

3.14 X 10-36 cm^l 

\lGeV J 



10^ < < 10^ GeV 



10« < < 10^' GeV 



(66) 



cr, 



CC 



0.484 



6.55 X 10-37 cm^( J \ 
5.20 X 10-36 cm'^ L^ly ) 



10^ <Etj< 10^ Ge\/ 



10« <Etj< 10^' GeV 



(67) 



0" 



VN 
NC 



3.04 X 10-37 



3.14 X 10-36 cm^ 



E^ 

1 GeV 
Etj \°-349 

IGel/y 



10^ < < 10« GeV 



10« < < 10^^ Gel^ 



(68) 



Our expressions are also to be compared with the NLO expressions for 
cross-sections of Ref. |2], calculated on the basis of CTEQ3-DIS parton 
distributions [83] . These cross-sections are well represented by simple power- 
law form for neutrino energies in the range 10^ GeV < Ei, < 10^^ GeV as 
given below: 



= 2.69 X lO-^'' cm' 



-36 „2 



Eu \ 

GeV J 



0.402 



(69) 



a^/c = 1.06 X lO-^'' cm' 



-36 „2 



1 Gev) 



(70) 
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rp_ X 0.404 

all = 2.53 X 10-^^ cm' [-^) (71) 

. 0.410 



a];- = 0.98 X 10-^^ cm' (^^^ j (72) 

Next, we compare our results with the NLO resuhs of Ref. [57], calculated 
on the basis of CTEQ4-DIS parton distributions [84J. For neutrino energies 
in the range 10^ GeV < < 10^^ GeV , the CTEQ4-DIS cross sections are 
given within 10% by 



P X 0.363 

afa = 5.53 X 10-^^ cm' (^^) (73) 



p s 0.363 

afa = 2.31 x 10^^^ cm' (74) 



P s 0.363 

= 7.84 X 10-^^ cm' ( ^) (75) 



p_ V 0.363 

4^ = 5.52 X 10-^^ cm' (76) 

^_ X 0.363 



<^ = 2.29 X 10-^^ cm' (^^j (77) 



all = 7.80 X 10-^« cm' (^) (78) 
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It is pertinent to note that CTEQ3-DIS parton distributions are more 
singular than CTEQ4-DIS parton distributions as x becomes ultra small i.e. 
as a; — )■ [SZ]. Around a; — )■ 0, CTEQ3 sea-quark distributions are presented 
by 

xq™\x) oc (79) 

On the other hand, the behaviour of CTEQ4 sea-quark distributions is 
well-presented by 

xq™\x) oc (80) 

Next we compare our analytical and numerical results with the analytical 
asymptotic approximation technique of Fiore et al as in Ref. [5] and their 
numerically obtained results [99] . 

We also compare our analytical and numerical results with the numeri- 
cal results of Rahul Basu et al \i\ obtained by using GRV 98 NLO parton 
distributions. 

We next find the ratio of o"^^ and and the ratio comes closer to the 
value 1 at the higher end of the neutrino (antineutrino) energy spectrum. 
This is because at higher energies (consequently at lower x), the contribution 
due to xFs part (valence part) becomes neghble in comparison to F2 part (sea 
and gluon part). 
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Energy cum vN interaction cross-section 




a^c icm ) 


^NC icm ) 


^tot icm ) 


10 


4.36 X 10-2^ 


7.24 X 10-39 


5.08 X 10-38 


102 


6.22 X 10-37 


1.08 X 10-37 


7.30 X 10-37 


10^ 


6.82 X 10-36 


1.22 X 10-36 


8.04 X 10-36 


10^ 


5.17 X 10-35 


1.00 X 10-35 


6.17 X 10-35 


10^ 


1.66 X 10-34 


3.59 X 10-35 


2.02 X 10-34 


10^ 


2.97 X 10-34 


6.82 X 10-35 


3.65 X 10-34 


10^ 


4.50 X 10-34 


1.05 X 10-34 


5.55 X 10-34 


10« 


6.68 X 10-34 


1.60 X 10-34 


8.28 X 10-34 


109 


1 X 10-33 


2.43 X 10-34 


1.24 X 10-33 


lO^o 


1.51 X 10-33 


3.70 X 10-34 


1.88 X 10-33 


10" 


2.31 X 10-33 


5.66 X 10-34 


2.88 X 10-33 


10^2 


3.54 X 10-33 


8.71 X 10-34 


4.41 X 10-33 



Table 1: Analytical estimation of neutrino-induced charged current, neutral 
current and total cross-sections in cm^ on the basis of eq. ( H9|) and eq. ( l50l) . 
considering the simple ansatz for the MRST f4 2004 input parton distribu- 
tions [5HI55] 
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Energy cum vN interaction cross-section 


TP (ncAi\ 

rLiu ) 


a^c icm ) 


^NC icm ) 


^tot icm ) 


10 


1.00 X 10-2^ 


2.46 X 10-39 


1.25 X 10-38 


102 


2.84 X 10-37 


6.01 X 10-38 


3.44 X 10-37 


10^ 


3.90 X 10-36 


8 X 10-37 


4.70 X 10-36 


10^ 


3.51 X 10-35 


7.40 X 10-36 


4.25 X 10-35 


10^ 


1.32 X 10-34 


3.01 X 10-35 


1.62 X 10-34 


10^ 


2.69 X 10-34 


6.30 X 10-35 


3.32 X 10-34 


10^ 


4.35 X 10-34 


1.02 X 10-34 


5.37 X 10-34 


10« 


6.61 X 10-34 


1.60 X 10-34 


8.21 X 10-34 


109 


9.96 X 10-34 


2.43 X 10-34 


1.24 X 10-33 


lO^o 


1.51 X 10-33 


3.70 X 10-34 


1.88 X 10-33 


10" 


2.31 X 10-33 


5.66 X 10-34 


2.88 X 10-33 


10^2 


3.54 X 10-33 


8.71 X 10-34 


4.41 X 10-33 



Table 2: Analytical estimation of antineutrino-induced charged current, neu- 
tral current and total cross-sections in cm^ on the basis of eq. ( ISTl) and 
eq. ( 152|) . considering the simple ansatz for the MRST f4 2004 input parton 
distributions [511 ES] 
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Energy cum vN interaction cross-section 




a^c icm ) 


^NC icm ) 


^tot icm ) 


10 


6.83 X 10-^^ 


1.40 X 10-38 


8.23 X 10-38 


102 


6.52 X 10-37 


1.32 X 10-37 


7.84 X 10-37 


10^ 


5.74 X 10-36 


1.16 X 10-36 


6.90 X 10-36 


10^ 


3.91 X 10-35 


8.10 X 10-36 


4.72 X 10-35 


10^ 


1.52 X 10-34 


3.16 X 10-35 


1.84 X 10-34 


10^ 


4.21 X 10-34 


8.53 X 10-35 


5.06 X 10-34 


10^ 


1.13 X 10-33 


2.21 X 10-34 


1.35 X 10-33 


10« 


3.04 X 10-33 


6.06 X 10-34 


3.65 X 10-33 


109 


8.07 X 10-33 


1.63 X 10-33 


9.70 X 10-33 


lO^o 


2.05 X 10-32 


4.17 X 10-33 


2.47 X 10-32 


10" 


4.95 X 10-32 


1.02 X 10-32 


5.97 X 10-32 


10^2 


1.14 X 10-31 


2.35 X 10-32 


1.38 X 10-31 



Table 3: Numerical estimation of neutrino-induced charged current, neutral 
current and total cross-sections in cm2 on the basis of eq. ( 16T1) and eq. ( l62l) . 
considering the GRV 94 parton distributions (which are dependent on both 
X and (52) [56j 
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Energy cum vN interaction cross-section 


TP (ncAi\ 

rLiu ) 


a^c icm ) 


^NC icm ) 


^tot icm ) 


10 


2.97 X 10-2^ 


2.98 X 10-39 


3.27 X 10-38 


102 


3.16 X 10-37 


3.61 X 10-38 


3.52 X 10-37 


10^ 


2.96 X 10-36 


3.56 X 10-37 


3.32 X 10-36 


10^ 


2.37 X 10-35 


3.30 X 10-36 


2.70 X 10-35 


10^ 


1.20 X 10-34 


2.06 X 10-35 


1.41 X 10-34 


10^ 


3.95 X 10-34 


7.54 X 10-35 


4.70 X 10-34 


10^ 


1.11 X 10-33 


2.21 X 10-34 


1.33 X 10-33 


10« 


3.04 X 10-33 


6.06 X 10-34 


3.64 X 10-33 


109 


8.07 X 10-33 


1.63 X 10-33 


9.70 X 10-33 


lO^o 


2.05 X 10-32 


4.17 X 10-33 


2.47 X 10-32 


10" 


4.95 X 10-32 


1.02 X 10-32 


5.97 X 10-32 


10^2 


1.14 X 10-31 


2.35 X 10-32 


1.38 X 10-31 



Table 4: Numerical estimation of antineutrino induced charged current, neu- 
tral current and total cross-sections in cm2 on the basis of eq. ( l63l) and eq. 
( I64l) . considering the GRV 94 parton distributions (which are dependent on 
both X and Q"^) (56] 
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Value of s cum ^^^^^ interaction cross-section (cm^) 


s (GeV2) 


a^^, Ref. [58] 


a^^, Ref. pj 


(T^c-; analytical 


(j'cc-i numerical 


10^ 


9.75 X 10-35 


1.03 X 10-34 


1.13 X 10-34 


9.25 X 10-35 


10^ 


4.13 X 10-34 


3.65 X 10-34 


2.41 X 10-34 


2.99 X 10-34 


10^ 


1.34 X 10-33 


1.11 X 10-33 


3.89 X 10-34 


8.26 X 10-34 


10^ 


3.75 X 10-33 


3.20 X 10-33 


5.87 X 10-34 


2.25 X 10-33 


109 


9.67 X 10-33 


8.61 X 10-33 


8.82 X 10-34 


6.04 X 10-33 


10^° 


2.33 X 10-32 


2.16 X 10-32 


1.33 X 10-33 


1.56 X 10-32 


10" 


5.34 X 10-32 


5.11 X 10-32 


2.03 X 10-33 


3.81 X 10-32 


10^2 


1.17 X 10-31 


1.14 X 10-31 


3.11 X 10-33 


8.89 X 10-32 


10^3 


2.45 X 10-31 


2.44 X 10-31 


4.78 X 10-33 


1.98 X 10-31 


10^^ 


4.97 X 10-31 


5.02 X 10-31 


7.36 X 10-33 


4.23 X 10-31 



Table 5: Comparison between average neutrino(antineutrino)-nucleon 
charged current cross section in cm2 according to the numerical results ob- 
tained in Ref. [99] (2nd column), the asymptotic approximation method of 
Ref. (3rd column), our analytical evaluation (4th column) obtained us- 
ing eq. fl49p as well as eq. (IHTj) and our numerical evaluation (5th column) 
obtained using eq. (|6T|) as well as eq. ( 163|1 . 
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Upper bounds on cross-section 


(GeV) 


Rice llOUUlOlj 


Agasa |102j 


J^'KRi' [103j 


FJ |104j 


FKRi' [103j 


FJ |104j 




1.2 X 10^3° 


2.8 X 10-31 


1.3 X 10-28 


2.4 X 10-29 


1010.5 


3.6 X 10-3° 


7.2 X 10-31 


No result. 


1.4 X 10-28 




1.9 X 10-29 


3.8 X 10-30 


No result. 


No result. 



Table 6: Results of the upper bounds on cross-section cm2 in Rice |100l 1101] 
and Agasa |102] experiment as per the estimation of neutrino flux of Ref. 
[103] and Ref. [TOl] . 



Our results versus Upper bounds on cross-section 


(GeV) 


Analytical 


Numerical 


Ref [2] 


Ref jl05] 


lO^o 


1.88 X 10-33 


2.47 X 10-32 


4.35 X 10-32 


2.24 X 10-32 


IQio.s 


2.32 X 10-33 


3.28 X 10-32 


No result. 


No result. 


lO^i 


2.88 X 10-33 


5.97 X 10-32 


1.02 X 10-31 


3.66 X 10-32 



Table 7: Our analytical and numerical results alongwith the results of the 
Ref. |2] and Ref. |105] to be compared with the experimental upper bounds 
on cross-section in cm2 as given in Table O 
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Sensitivity of PAO 


(GeV) 


FKRT |lU3j 


PJ [TM] 




1.0 X 10^30 _ ^ g X 10-30 


3.1 X 10-31 - 5.6 X 10-31 




4.2 X 10-30 - 8.7 X 10-30 


1.1 X 10-30 - 2.1 X 10-30 


10" 


2.7 X 10-29 - 7.8 X 10-29 


6.2 X 10-30 - 1.5 X 10-29 



Table 8: Sensitivity of Pierre Auger Observatory (PAO) at 95% confidence 
limit in cm^. This is derived with the assumption that no deeply developing 
showers above Standard Model background has been observed in 5 years of 
operation. This can be compared with the results given in Table [3 
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Energy cum total neutrino-nucleon interaction cross-section 


(GeV2) 


Ref. P 


Ref. [57] 


<otah analytical 


(^totah numerical 


10^ 


2.48 X 10-33 


2.72 X 10-33 


5.55 X 10-34 


1.35 X 10-33 


108 


5.59 X 10^33 


6.29 X 10-33 


8.28 X 10-34 


3.65 X 10-33 


10^ 


1.26 X 10-32 


1.45 X 10-32 


1.24 X 10-33 


9.70 X 10-33 


10^0 


2.83 X 10-32 


3.34 X 10-32 


1.88 X 10-33 


2.47 X 10-32 


10" 


6.37 X 10-32 


7.71 X 10-32 


2.88 X 10-33 


5.97 X 10-32 



Table 9: Comparison between the neutrino-nucleon total cross section in 
cm2 according to the numerical results obtained in Ref. [Tj (2nd column), 
the numerical method of Ref. [57] (3rd column), our analytical evaluation 
(4th column) obtained using eq. (H9|) as well as eq. ( 150|) and our numerical 
evaluation (5th column) obtained using eq. as well as eq. (1^ . 
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Energy versus — ^ 




E.iv) (GeV^) 


analytical 


numerical 


Ref |106j 


10 


4.36 


2.30 


NO. 


102 


2.19 


2.06 


z.Ud 


103 


1.75 


1.94 


1 so 


10^ 


1.47 


1.65 


1 c;o 
l.dZ 


10^ 


1.26 


1.27 


1 1 o 
i.iy 


10^ 


1.10 


1.07 


i.Uo 


10^ 


1.03 


1.02 


i\o. 


108 


1.01 


1.00 


No. 


109 


1.00 


1.00 


No. 


10^° 


1.00 


1.00 


No. 


10" 


1.00 


1.00 


No. 


10^2 


1.00 


1.00 


No. 



Table 10: The ratio of neutrino-nucleon charged current cross-section and 
antineutrino-nucleon charged current cross-section evaluated by our analyti- 
cal and numerical method is plotted as a function of neutrino (antineutrino) 
energy {E^(j^^\. This is compared with the ratio obtained from Ref [106] 
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In Figure El we plot our analytical LO results (Magenta coloured dash 
dot line) for g^qq [cvn^^ eq. (l49l) (Table [1]). We then plot exact numerical LO 
results (Blue coloured dash line) calculated by us eq. fl^T]) (Table [3]). In the 
same figure, we also plot NLO results of Ref. eq. f l^ (Pink coloured dot 
line) and that of Ref. [57] (Green coloured solid line) eq. fl69|l as well as Ref. 
[2] eq. (175]) (Wine coloured dash dot dot line). 

In Figure [TJ we compare our LO predictions (Magenta coloured dash dot 
line) for g^j^q [cvp?^ eq. (l50|l (Tabled]), with the exact numerical LO results 
calculated by us eq. fl52]) (Blue coloured dash line) (Table E]). The NLO 
results of Ref. |3] eq. fl66]) (Pink coloured dot line) and that of Ref. [57] 
(Green coloured solid line) eq. (17D]) as well as Ref. [2] eq. (|71]) (Wine coloured 
dash dot dot line) are then plotted to have a relative comparison. 

Similarly in Figure El we compare our LO results for g^qq [cvn}^ (Magenta 
coloured dash dot line) eq. ( ISTl) (Table [2]) with the exact numerical LO 
results calculated by us eq. f l63|) (Blue coloured dash line) (Table H]), the 
NLO results of Ref. [3] eq. ([67D (Pink coloured dot line), Ref. |7] (Green 
coloured solid line) eq. ( I7T]) as well as Ref. [2] eq. ( 17^ (Wine coloured dash 
dot dot line). 

In Figure [9l we plot our LO results for cr)^^ (cm^) (Magenta coloured dash 
dot line) eq. ( I52l) (Table [2]) and then compare with the exact numerical LO 
results calculated by us eq. flMj) (Blue coloured dash line) (Table H]). For 
relative comparison, the NLO results of Ref. [3] eq. (l68ll (Pink coloured dot 
line), Ref. [57] (Green coloured solid line) eq. ( rf2l) as well as Ref. p] eq. 



45 



([77]) (Wine coloured dash dot dot line) are also plotted. 

In Figure [TOl we compare the value of obtained by using our ana- 
lytical (Magenta coloured dash dot line) (using eq. ( H9l) as well as eq. ( l5Ti) ) 
and numerical techniques (Wine coloured dot line) (using eq. ( I6T1) as well 
as eq. fl63|) ) with the numerical results obtained in Ref. [99] (Blue coloured 
solid line), the asymptotic approximation method of Ref. ^ (Red coloured 
dash line) (Table |5]). It is pertinent to note that the factor is defined as: 

uN I VN 
-=vN _ ^CC + "cc 

^cc - ^ • 

In Figure [m we compare between the neutrino- nucleon total cross section 

in cm^ according to the numerical results obtained in Ref. [1] (Blue coloured 
solid line), the numerical method of Ref. [57] (Magenta coloured dash line), 
our analytical evaluation (Olive coloured dot line) and our numerical evalu- 
ation based on GRV94 distribution (Red coloured dash dot dot line) (Table 



In Figure [121 the ratios of neutrino-nucleon charged current cross-section 
and antineutrino-nucleon charged current cross-section for both of 



our analytical (Blue coloured solid line) and numerical methods (Red coloured 
dash line) respectively are plotted as a function of neutrino (antineutrino) 



energy [E^i^y^. The ratio obtained from Ref [106] (Magenta coloured dash 
dot line) is compared with our analytical and numerical results (Table fTOj) . 



[9]). 
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3.1 Comparative study and limitations: 

Let us make a comparative study of the Figures (PfTTl). As expected and 
clearly demonstrated in these figures, our analytical result indicates a sharp 
rise for < 10^ GeV. This is due to the absence of propagator effect in low 
energy regime [69]. For E^, > 10^ GeV, the dampening due to the propagator 
effect takes place. Here we find our numerically determined result in sufficient 
agreement with the results obtained by the above-mentioned authors PQ El 
[31 [571 [29] • As far as our analytical result is concerned, we find it to match 
better with other results at lower end of the energy spectrum (10^ < E^(p) < 
10^) of the UHE neutrino rather than the higher end (10^ < -^1/(77) < 10^^). 
At the higher end, we observe dampening of the cross-section than what is 
expected from the other numerical results etc. mentioned beforehand. The 
quantitative difference at the higher end is presumably due to the following 
factors: 

(1) We have not extrapolated MRST distribution for values of x below 
10^^, assuming that our approximate analytical solutions of DGLAP equa- 
tions for singlet structure function F2{x,t) as well as non-singlet structure 
function xF^ remains valid in the UHE regime. 

(2) We have considered solutions at tree level. Replacing LO distributions 
by NLO distributions and with the appropriate addition of convolutions with 
the fermionic Wilson coefficient Cg and the NLO contribution Cg ® g [93] , 
we may expect somewhat better results. 

(3) We have neglected the finite x corrections coming from the higher 
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dF^^(x,t) 

derivatives of — ^ — in Taylor approximation of DGLAP equations in 

ox 

case of t-evolution of xF^ at one loop level [52]. Similarly, we have neglected 

dF^ix, t) 

the finite x corrections coming from the higher derivatives of — ^ — and 

ox 

dG(x,t) 

— ^ — in Taylor approximation of DGLAP equations in case of t-evolution 

Ox 

of F2 at one loop level [53] . 

(4) We have considered the contribution due to light quarks only and ne- 
glected the contribution due to heavy quarks. For Ei, > 10^ GeV, preferably 
for El, — 7- 10^^ GeV, NLO massive tb contribution will introduce an 0{as) 
correction of about 20% to the LO contribution of cr^^ . 

We find that the blue coloured dashed line drawn on the basis of our 
numerical method has got the almost same slope as that of the other lines 
drawn on the basis other numerical methods as seen in Figures ([Ml]), proving 
the excellence of our numerical results. A little lateral shift is presumably 
because of the following reasons: 

(1) We have done calculations at leading order involving LO parton dis- 
tributions, whereas the calculations done by other authors at next-to-leading 
order are based on NLO parton distributions. Obviously a better result is 
expected at NNLO level. 

(2) It has been found that there had been a 10% discrepancies between 
GRV 94 distribution [56] adopted by us for numerical calculation and the 
precision measurements at HERA done at a later stage (in 1996 and 1997). 
This implies that GRV 94 distribution has restricted accuracy margin. So 
Gluck et al realised the necessity of the fine tuning of the input parameter n 
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and/or f{x, /i^) both at LO and NLO level which they did and removed these 
10% discrepencies by providing the GRV 98 parton distributions [63j- It is 
pertinent to remind that NLO calculations of Gluck et al were based mainly 
on GRV 98 parton distributions, which was almost free from any uncertainty 
of radiative predictions because of timely fine tuning. For example, they 
fine-tuned /i^q from the value 0.23 GeV^ to 0.26 GeV^ at one loop level. 
Similarly there was a fine tuning at two loop level from from the value 0.34 
GeV2 to 0.40 GeV2 [561 [63] • 

(3) GRV 94 predictions adopted by us are valid in the small x region 
10^^ < a; < 1. We did our numerical calculations assuming the validity of 
these distributions in the UHE regime i.e. 10^^^ < a; < 10^^, which is obvi- 
ously not true. Uncertainty and error are definite to crop up (however small 
it may be) due to the this non-extrapolation. It is reminded that Gluck et al 
later extended their perturbatively stable parameterfree dynamical predic- 
tions to the extremely small-x region 10"^ < x < 10~^ [63] . 

(4) We have not taken into consideration the contribution due to heavy 
quarks. 

3.2 Upper bounds on the neutrino-nucleon inelastic 
cross-sections: 

Model independent upper bounds on the neutrino-nucleon inelastic cross- 
sections derived from the RICE [TOOl [M] and AGASA ^0^ collaboration 
search results by exploiting the cosmogenic neutrino flux estimates of Ref. 
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|1U3] and Ref. |1U4] are already known. This is shown in Table E] |107] . 
Such an upper bound arises due to non-observation of events triggered off 
by UHE neutrinos as reported by RICE and AGASA. We compare our an- 
alytical as well as numerically determined results for the Ultra High Energy 
neutrino-nucleon cross-section in leading order with the above-mentioned up- 
per bounds at three different neutrino energies E^^ = 10^^, 10^*^'^, 10^^ GeV. 
We find that our results do not violate the experimental upper bounds im- 
posed on the neutrino-nucleon interaction cross-section and are well within it 
as clearly seen in Table Eland Table [71 In Table |H1 Sensitivity of Pierre Auger 
Observatory (PAO) at 95% confidence limit in cm^ is shown |107j . This has 
been derived with the assumption that no deeply developing showers above 
Standard Model background has been observed in 5 years of operation. This 
can be compared with the results given in Table [71 

In Figure [131 we plot our analytical result (Wine coloured dot line) 
and numerical (exact) result (Magenta coloured dash dot line) for neutrino- 
nucleon cross-section derived from Standard Model. This is compared with 
the upper bounds on the neutrino-nucleon cross-sections derived from the 
RICE collaboration search results by exploiting the cosmogenic neutrino flux 
estimates of FKRT |103] (Red coloured Solid line) and PJ |104] (blue coloured 
dash line). 
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3.2.1 Beyond the standard model scenario: 

As per Beyond the Standard model scenario, neutrinos with energies around 
lO^*^'^ GeV may behave as strongly interacting as protons, which has the ca- 
pacity of inducing vertical air showers at high altitudes in the atmosphere. 
Since there is a good match of cosmogenic neutrino flux with the UHE cosmic 
ray spectrum above Greison-Zatsepin-Kuzmin (GZK) energy {Egzk ~ 10^°'^ 
GeV), so one can appreciate this 'strongly interacting' characteristics of cos- 
mogenic neutrinos which can solve the GZK puzzle. If the strongly inter- 
acting neutrino scenario is considered for Ey >z 10^^ GeV, which leads to 
the generation of vertical showers, then non-observation of quasi-horizontal 
air showers by AGASA would not be able to constrain the neutrino-nucleon 
cross-section the way it did. 

4 Comments and conclusions: 

To conclude, we have shown that our analytical evaluation of neutrino (antineutrino)- 
nucleon cross-section tallies well in selected part of of the neutrino(antineutrino) 
energy spectrum. Particularly it tallies well with the result obtained at the 
lower end of the energy spectrum, but at the higher end there is dampening 
much more than what is expected due to propagator effect in case of nu- 
merical solutions by other authors [11 |2l |3l |5l [571 [99] • The reason for such 
differences at higher end of the neutrino(antineutrino) energy spectrum has 
been explained in subsection 13.11 On the other hand, our LO numerical 
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result tallies quite well with the NLO results of other authors both at the 
lower as well as higher ends of the neutrino(antineutrino) energy spectrum. 
The possible reasons for minute lateral shift has been explained adequetly in 
subsection 13.11 
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Figure 6: UHE neutrino- nucleoli charged current cross-section cr^^ (cm^) 
obtained from our analytic LO expression eq. ( l49i) (Magenta coloured dash 
dot line). This is compared with our numerical LO expression eq. ( 16T1) (Blue 
coloured dash line). In the same figure, we also plot NLO results of Ref. [3] 
eq. (l65l) (Pink coloured dot line) and that of Ref. [57j (Green coloured solid 
line) eq. (169|1 as well as Ref. [2] eq. (173|l (Wine coloured dash dot dot line) 
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Figure 7: UHE neutrino-nucleon neutral current cross-section a'^ (cm^) ob- 
tained from our analytic LO expression eq. f lSU]) (Magenta coloured dash dot 
line). This is compared with our numerical LO expression eq. f E2]) (Blue 
coloured dash line). In the same figure, we also plot NLO results of Ref. [3] 
eq. (!66|) (Pink coloured dot line) and that of Ref. [57j (Green coloured solid 
line) eq. (ITOjl as well as Ref. [2] eq. (17i|l (Wine coloured dash dot dot line) 
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Figure 8: UHE anti-neutrino nucleon charged current cross-section {cm?) 
obtained from our analytic LO expression eq. flSTl) (Magenta coloured dash 
dot line). This is compared with our numerical LO expression eq. (jHS]) (Blue 
coloured dash line). In the same figure, we also plot NLO results of Ref. [3] 
eq. (1671) (Pink coloured dot line) and that of Ref. [57] (Green coloured solid 
line) eq. (17T]) as well as Ref. [2] eq. (176|) (Wine coloured dash dot dot line). 
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Figure 9: UHE anti-neutrino nucleon neutral current cross-section o^j^q {cm'^) 
obtained from our analytic LO expression eq. f l52|) (Magenta coloured dash 
dot line). This is compared with our numerical LO expression eq. f lM|) (Blue 
coloured dash line). In the same figure, we also plot NLO results of Ref. [3] 
eq. (168!) (Pink coloured dot line) and that of Ref. [57] (Green coloured solid 
line) eq. (!72|) as well as Ref. [2] eq. (!77|) (Wine coloured dash dot dot line). 
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Figure 10: We compare between average neutrino(antineutrino)-nucleon 
charged current cross section in cm^ according to the numerical results ob- 
tained in Ref. [HS] (Blue coloured solid line), the asymptotic approximation 
method of Ref. [3] (Red coloured dash line), our analytical evaluation (Ma- 
genta coloured dash dot line) obtained using eq. (H9|l as well as eq. (!5T|) and 
our numerical evaluation obtained (Wine coloured dot hne) using eq. ( I6T!) 
as well as eq. ( 163|) 
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Figure 11: We plot the neutrino-nucleon total cross sections in cm^ versus 
neutrino energy to compare the numerical results obtained in Ref. [1] (Blue 
coloured solid line), the numerical method of Ref. [57] (Magenta coloured 
dash hne), our analytical evaluation obtained using eq. fH^ (Olive coloured 
dot line) as well as eq. ( 150|1 and our numerical evaluation obtained using eq. 
(!6T|) (Red coloured dash dot dot line) as well as eq. (!62|) 
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Figure 12: We plot the ratio of neutrino- nucleon charged current cross-section 

f a"^ \ 

and antineutrino-nucleon charged current cross-section of our analyt- 

X'^cc J 

ical (Table 6.1 and 6.2) (Blue coloured solid line) and numerical method 
(Table 6.3 and 6.4) (Red coloured dash line). The ratio obtained from Ref 
[106] (Magenta coloured dash dot line) versus neutrino (antineutrino) energy 
{E^ijj'^ is also plotted for comparison. 
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Figure 13: Model independent upper bounds on the neutrino-nucleon inelas- 
tic total cross-sections derived from the RICE |100l 1101] collaboration search 
results by exploiting the cosmogenic neutrino flux estimates of Ref. |103] 
and Ref. |104j are plotted against neutrino energy. Our analytical as well as 
numerically (exact) determined results from Standard Model for the Ultra 
High Energy neutrino-nucleon total cross-section in leading order are com- 
pared with the above-mentioned upper bounds at three different neutrino 
energies = 10^°, 10l°•^ lO^^ GeV. 
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